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Abstract 

We  present  autoregressive  (AR)  and  autoregressive  moving  average  (ARMA) 

I 

processes  with  bivariate  exponential  (BE)  and  bivariate  geometric  (BG)  distri¬ 
butions.  The  theory  of  positive  dependence  is  used  to  show  that  in  various 
cases,  the  BEAR,  BGAR,  BEARMA,  and  BGARMA  models  consist  of  associated  random 
variables.  We  discuss  special  cases  of  the  BEAR  and  BGAR  processes  in  which  | 

I 

the  bivariate  processes  are  stationary  and  have  well  known  bivariate  exponen-  j 

! 

tial  and  geometric  distributions.  | 

i 


1.  Introduction  and  Sunmiary 

A  primary  stationary  model  in  time  series  analysis  is  the  pxl  linear  pro¬ 
cess  given  by: 

CO 

(1.1)  X(n)  -  I  A(j)e(n-j),  n  -  0,±1.±2 . 

j— «0 

where  A( j ) ,  j  -  0,±1,±2 .  is  a  sequence  of  pxp  parameter  matrices  such  that 

CO 

X  l|A(j)||  <  ®,  and  £(n) ,  n  -  0,±1,±2,...,  is  a  sequence  of  uncorrelated 
J— « 

pxl  random  vectors  with  mean  zero  and  common  covariance  matrix.  It  is  well 
known  that  (1.1)  includes  the  stationary  vector  autoregressive  (AR)  process 
and  the  stationary  and  invertible  vector  autoregressive  moving  average  (AEMA) 
process.  However,  in  some  physical  situations  where  the  random  vectors  X(n) 
are  either  positive  or  discrete,  the  preceeding  assumptions  on  the  £(n) 
sequence  are  inappropriate  (see  Lewis  (1980),  p.  152). 

Several  researchers,  addressing  themselves  to  this  problem,  have  con¬ 
structed  univariate  stationary  AR  type  models  and  stationary  ARMA  type  models 
where  the  random  variables  X(n)  have  exponential  or  gamma  distributions,  and 
discrete  models  where  X(n)  assumes  values  in  a  common  set.  Lawrence  and  Lewis 
(1977,  1980,  1981)  and  Jacobs  and  Lewis  (1977)  present  stationary  AR  and  ARMA 
type  models  where  the  random  variables  X(n)  have  exponential  distributions; 
Gaver  and  Lewis  (1980)  consider  stationary  ARMA  type  models  where  the  random 
variables  X(n)  have  gamma  distributions.  Jacobs  and  Lewis  (1978a, b,  1983) 
construct  ARMA  type  models  where  the  random  variables  X(n)  are  discrete  and 
assvime  values  in  a  common  finite  set.  The  aforementioned  models  have  been 


used  in  the  various  fields  of  applied  probability  and  time  series  analysis, 
for  example,  these  models  have  been  used  to  model  and  analyze  univariate  point 
processes  with  correlated  service  and  correlated  interarrival  times  (see 


Jacobs  (1978)).  Details  concerning  bivariate  exponential  and  geometric  MA 
type  processes  and  the  corresponding  point  processes  may  be  found  in  Langberg 
and  Stoffer  (1985). 

In  this  paper  we  present  two  classes  of  AR  and  of  ARMA  type  sequences  of 
bivariate  random  vectors.  The  first  class  has  exponential  marginals  while  the 
sec'  .d  class  has  geometric  marginals.  We  denote  the  first  [second]  class  of 
models  as  BEAR(m)  [BGAR(m)]  and  BEARMA(raj^ ,m2)  [BGARMA(mj^  ,m2)  ]  for  bivariate 
exponential  [geometric]  autoregressive,  order  m,  and  autoregressive  moving 
average,  order  (m^,m2),  respectively,  where  m  and  (mj^,m2)  parametrize  the 
order  of  the  dependence  on  the  past.  We  use  the  theory  of  positive  dependence 
to  show  that  in  a  variety  of  cases  the  classes  of  sequences  are  associated. 

In  Section  2  we  define  the  bivariate  exponential  and  geometric  distribu¬ 
tions  which  are  the  underlying  distributions  of  our  two  classes,  and  present  a 
variety  of  examples  of  such  distributions.  Furthermore,  in  Section  2  we 
define  the  concept  of  association  and  present  a  variety  of  bivariate  exponen¬ 
tial  and  geometric  distributions  that  are  associated.  We  conclude  Section  2 
by  describing  the  bivariate  dependence  mechanisms  which  are  used  in  generating 
the  various  models.  In  Section  3  we  construct  the  general  BEAR(m)  and 
BGAR(m)  models  showing  that  the  sequences  have  bivariate  exponential  and 
bivariate  geometric  distributions,  respectively.  Also,  we  discuss  the  auto¬ 
correlation  structure  of  the  two  classes  of  sequences.  In  Section  4  we  con¬ 
sider  special  cases  of  the  BEAR(l)  and  BGAR(l)  sequences.  We  show  that 
defined  appropriately,  the  bivariate  processes  are  stationary,  and  obtain  well 
known  bivariate  exponential  and  geometric  distributions.  Finally,  in  Section 
5,  we  present  the  BEARMA  (mj^,m2)  and  BGARMA(mj^ ,m2)  models.  We  conclude  by 
describing  the  association  properties  of  the  sequences  and  discussing  how  to 
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utilize  association  to  obtain  some  probability  bounds  and  moment  inequalities 


for  the  bivariate  processes  and  the  corresponding  point  processes. 


2.  Preliminaries 


In  this  section  we  present  definitions  and  prove  some  basic  results  to  be 


used  in  the  sequel.  First,  we  give  definitions  of  bivariate  geometric  and 


exponential  distributions  and  provide  some  examples.  Then  we  present  the  con¬ 


cept  of  association  and  give  some  examples.  Finally,  we  discuss  some  bivari¬ 


ate  dependence  mechanisms . 


First,  we  present  a  definition  of  a  bivariate  geometric  distribution. 


Definition  2.1.  Let  M,N  be  random  variables  assuming  values  in  the  set 


(1,2,...).  We  say  that  (M,N)  has  a  bivariate  geometric  distribution  if  M  and 


N  have  geometric  distributions. 


Examples  2.2.  (a)  Let  N  be  geometric.  Then  (N,N)  is  bivariate  geometric. 


(b)  Let  M  and  N  be  Independent  geometric  random  variables,  then  (M,N)  is 


bivariate  geometric.  (c)  Let  Nj^,N2,N2  be  independent  geometric  random  vari¬ 
ables,  and  put  M  -  {min(Nj^  ,N2)  ) ,  N  -  (min(N2  .N^)  ) .  Then  (M,N)  has  the  Esary- 


Harshall  (1974)  bivariate  geometric  distribution.  (d)  Let  Poo’^Ol’^lO’^ll 


in  [0,1]  such  that  (i)  Poo+Poi'^PlO'^Pll  “  Pqi'^PiI  ^  PlO'^Pll  ^ 

and  let  M,N  be  random  variables  assuming  values  in  the  set  (1,2,...)  deter¬ 


mined  by; 


(2.3)  P(M  >  m,  N  >  n)  - 


™  r  ,n*in 

Pll[Poi+Pii]  -  n  >  m. 


Xi  f  .  1  ®  ^  ^  4  /\ 

PlllPlO+Piil  .  n  <  m.  m.n  -  1,2, 


‘.♦'I 


I 

>*l 
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Then  (M,N)  has  the  Block  (1977)  fundamental  bivariate  geometric  distribution 
(see  also  Block  and  Paulson  (1984)).  (e)  Let  (Mj^,M2j  be  bivariate  geometric 

and  let  (Nj^(j )  ,N2(j ) ) ,  j  -  1,2,...,  be  an  iid  sequence  of  random  vectors  with 
bivariate  geometric  distributions  which  are  independent  of  (Mj^,M2j.  Then 
Ml  M2 

(  Z  N..(j),  X  N-(j))  has  a  bivariate  geometric  distribution  (cf.  Lemma  2.14). 

j-1  j-1 

In  the  following  remark  we  show  how  some  of  the  bivariate  geometric  dis¬ 
tributions  are  particular  cases  of  Example  2. 2d.  Other  examples  are  given  in 
Remark  4.10. 

Remarks  2.4.  (a)  Let  P]^o”^Ol”^  equation  (2.3).  Then  we  obtain  the 

distribution  introduced  in  Example  2.2a.  (b)  Let  p^^^  -  ^^11^^01^ 

(2.3).  Then  we  obtain  the  bivariate  geometric  distribution  introduced  in 

Example  2.2b.  (c)  Let  Pii><Pii+Pio^ ^Pll'*’Poi^  N^,N2,Nj  be 

independent  geometric  random  variables  with  parameters  * 

^ll^Pll'^PlO^*^’  Pll^^ll'^’PlO^  ^Pll'‘’P01^’  'Respectively.  Put  M  -  {min(Nj^,N2) } 
and  N  -  (min(N2 ,N2) ) .  Then  (M,N)  is  stochastically  equal  to  the  Esary- 
Harshall  bivariate  geometric  distribution  given  in  Example  2.2c. 

Next,  we  present  a  definition  of  a  bivariate  exponential  distribution. 

Definition  2.5  Let  Ej^,E2  be  random  variables  assuming  values  in  (0,®).  We  say 
that  (E^,E2)  has  a  bivariate  exponential  distribution  if  E^  and  E2  have 
exponential  distributions. 


Examples  2.6  (a)  Let  E  be  exponential.  Then  (E,E)  is  bivariate  exponential, 
(b)  Let  Ej^,E2  be  independent  exponentials.  Then  (Ej^,E2)  has  a  bivariate 
exponential  distribution.  (c)  Let  Xj^,X2,Xj  be  independent  exponentials  and 
put  Ej^  -  {mln(X^ ,X2) ) ,  E2  -  {min(X2 .X^^ ) .  Then  (E^,E2)  has  the  Marshall-Olkin 


(1967)  bivariate  exponential  distribution.  (d)  Let  (M,N)  have  a  bivariate 


geometric  distribution  and  let  (E^(j ) ,E2(j ) ) ,  j  -  1,2 .  be  an  iid  sequence 

of  random  vectors  with  bivariate  exponential  distributions,  independent  of  M 


M  N 

and  N.  Then  (  J.  Z  a  bivariate  exponential  distribution 

j-1  j-1 


(cf.  Lemma  2.14).  (e)  Let  0  <  <»  <  1.  then  (E^,E2)  determined  by 

P{Ej^  >  X,  E2  >  y)  “  exp{ -x-y-axy)  ,  x,y  >  0,  has  a  Gumbel  (1960)  bivariate 

exponential  distribution.  (f)  Let  |a|  <  1.  Then  (Ej^,E2)  determined  by 

P{Ej^  <  X,  E2  <  y)  -  (1-e  ’^)(l-e  (1+ae  *  ^)  ,  x,y  >  0,  has  a  bivariate  Gumbel 

(1960)  exponential  distribution.  (g)  Let  a  >  1.  Then  (Ej^,E2)  determined  by 


.  a  a. 1/a 

P{E^  >  X,  E2  >  y)  -  e  ^  ^  ,  x,y  >  0,  is  bivariate  exponential.  (h)  Let 


(X,Y)  be  a  random  vector  with  continuous  marginal  distributions  F  and  G, 
respectively.  Then  the  random  vector  (-ln[l-F(X) ] ,  -ln[l-G(Y)))  is  bivariate 
exponential . 


Example  2.6(d)  has  been  used  by  several  researchers  to  generate  bivariate 
distributions  (for  example  Arnold  (1975),  Downton  (1970),  and  Hawkes  (1972)  to 
mention  a  few) .  In  the  following  remarks  we  illustrate  how  some  of  the 
bivariate  exponential  distributions  are  obtained  from  Example  2.6(d).  Other 
examples  are  given  in  Remark  4,5. 


Remarks  2.7.  (a)  M  -  N  and  let  Ej^(j),  E2(j)  be  independent  exponentials, 

j  -  1,2,....  Then  we  obtain  the  distribution  introduced  by  Downton  (1970). 

(b)  Let  (M,N)  be  as  in  Example  2.2(d)  and  let  E^(j),  E2(j)  be  independent 
exponentials,  j  -  1,2,....  Then  we  obtain  the  bivariate  exponential  distribu¬ 
tion  introduced  by  Hawkes  (1972)  and  Paulson  (1973).  (c)  Let  (M,N)  be  as  in 

Example  2.2(c)  and  let  Ej^(j)  -  E2(j),  j  -  1,2,,...  Then  we  obtain  the 
Marshall -Olkin  (1967)  distribution  given  in  Example  2.6(c)  (for  details  see 
Marshall -Olkin  (1967)). 


Next,  we  present  a  concept  of  positive  dependence. 

Definition  2.8  Let  T  -  • • • ’^n^’  ””  1.2,...,  be  a  multivariate  random  vec¬ 
tor.  We  say  that  the  random  variables  T, . T  are  associated  if 

cov(f (T) , g(T) )  >  0  for  all  f  and  g  monotonically  nondecreasing  in  each  argu¬ 
ment,  such  that  the  expectations  exist. 

Remarks  2.9.  (a)  Note  that  independent  random  variables  are  associated 

and  that  nondecreasing  functions  of  associated  random  variables  are  associated 
(cf.  Barlow  and  Proschan  (1975)  pp.  30-31).  Thus  the  components  of  the  vector 
given  in  Example  2.2(c)  and  the  components  of  the  vector  given  in  Example 
2.6(c)  are  associated.  (b)  Let  (Ej^,E2)  be  as  in  Example  2.6(e),  with  a  >  0, 
or  as  in  2.6(f)  with  -1  <  a  <  0.  Since  P{E^  >  x,  E2  >  y)  <  P{E^  >  x)P(E2  >  y) 
for  x,y  >  0,  E^  and  E2  are  not  associated.  (c)  Let  (X,Y)  be  as  in  Example 
2.6(h).  Then  -ln[l-F(X)]  and  -ln[l-G(Y)]  are  associated  if  and  only  if  X  and 
Y  are  associated  (cf.  Barlow  and  Proschan  (1975),  Proposition  3,  p.  30).  (d) 

The  components  of  the  bivariate  geometric  distribution  given  in  Example  2.2(e) 
are  associated  provided  and  M2,  and  Nj^(l)  and  N2(l)  are  associated  (cf. 
Langberg  and  Stoffer  (1985),  Lemma  2.10).  (e)  The  components  of  the  bivariate 

exponential  distribution  given  in  Example  2.6(d)  are  associated  provided  that 
M  and  N,  and  E^(l)  and  E2(l)  are  associated  by  the  same  reasoning  as  in  Remark 
2.9(d). 

We  are  now  ready  to  discuss  the  various  dependence  mechanisms  used  in 
obtaining  bivariate  exponential  and  geometric  distributions.  It  turns  out 
that  many  of  these  mechanisms  are  related  and  we  describe  these  relationships . 

St 

The  notation  X  -  Y  will  mean  that  X  and  Y  are  random  variables  (or  vectors) 
with  the  same  distribution. 
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Lemma  2.W.  (Random  Mixing)  (a)  Let  (X^.X^)  and  be  independent  random 


St 


St 


vectors  with  exponential  marginals  where  X^  -  Z^,  X2  -  Z2 ,  and  (X^,X2), 
(Z^,Z2)  have  mean  vector  (Aj^^,A2^).  Let  (I^,l2)  be  a  bivariate  Bernoulli  ran¬ 
dom  vector  independent  of  (X^,X2),  (Zj^,Z2).  Also,  assume 


(2.11)  P{Ij^  -  i.  I2  -  j)  -  P^j  i.j  -  0,1 


such  that  \  Pj^j  -  1;  1-w^  -  P^o‘''Pll  ^“^2  ~  ^ 

^  I  j 


dom  vector  given  by 


St 


(2.12)  ^2^2^  ^2^'^ 


has  a  bivariate  exponential  distribution  with  the  same  marginals  as  (X^,X2) 


vec- 


and  (Zj^,Z2).  (b)  Let  (M^,M2)  and  (K^,K2)  be  independent  geometric  random 

tors  such  that  for  k  -  1,2 .  P{Mj^  -  k)  -  (p^^/x^)  (l-p^^^w^^)^  ^  and  P{Kj^  -  k) 

-  P^^.p  )^'^.  i  “  1.2,  with  0  <  Pj^.P2  <  1-  bet  (1^,12)  be  defined  by  (2.11) 


and  be  independent  of  (M^,M2),  (Kj^,K2).  Then  a  random  vector  given  by 


St 


(2.13)  (G^,G2)  -  (I^K^.  I2K2)  +  (M^,M2) 


has  the  same  marginals  as  (K^,K2). 


Proof:  Parts  (a)  and  (b)  follow  easily  by  computing  the  marginal  characteris¬ 
tic  functions.  || 


Raftery  (1985)  gives  a  special  case  of  Lemma  2.10a  where  ■  Z2 . 


Lemma  2 . 14  (Random  Summation) .  Let  (Xj^j,X2j)  [(M^j,M2j)],  j  -  1,2 .  be  iid 


bivariate  exponential  [geometric]  random  vectors  with  mean  vector 

(ir^/A^,  W2/A2)  [(Wj^/Pj^,  jr2/P2)].  0  <  ’^^^.^2  <  1-  Let  (N^,N2)  have  a  bivariate 


),>jrKk’>2 


m 
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■1  -1. 


geometric  distribution  with  mean  vector  )  and  be  independent  of  the 

(X^j,X2j)  and  the  Then  random  vectors  given  by 


St 


(2.15)  (Yn.Y.)  -  (  I  X  I  X  ) 

j-1  J  J-1 


St  ^1 

(2.16)  -  (X  M  I  M  ) 

j-1  J-1 

have  bivariate  exponential  and  bivariate  geometric  distributions  with  mean 
vectors  (A^^,A2^)  and  (p^^,P2^).  respectively. 

Proof :  The  proof  follows  easily  by  computing  the  marginal  characteristic  func¬ 
tions  .  I  I 

Finally,  we  describe  how  random  mixing  and  random  summation  are  related 
in  the  following  lemma  and  corollary.  The  connection  between  the  two  concepts 
is  a  key  element  in  the  development  of  the  bivariate  exponential  and  geometric 
AR  and  ARMA  models  discussed  in  Sections  3,4,  and  5. 


Lemma  2.1J^.  Let  (X^j,X2j)  and  (Mj^j,M2j),  j  -  1,2 . be  as  defined  in  Lemma 

2.14.  Let  (N^,N2)  have  the  bivariate  geometric  distribution  given  by  (2.3) 
with  ”  P^o^^ll  ^  ”  ^01^^11  ^  independent  of  the 

(X^j,X2j)  and  (M^j,M2j).  Furthermore,  let  (Xj^,X2).  (Ij^,l2),  (Mj^,M2)  and 

Ni  N2 

(K..,K2)  be  as  defined  in  Lemma  2.10.  Then  (Y.,  ,Y_)  -  (  X  X.,  ,  X  X„  )  has  the 

h  ^2 

representation  (I^Z^,l2Z2)  +  (x^^X^ ,^2X2) .  and  (^^^,62)  -  (  X  ^  M2 ^ )  has 

the  representation  12^2^ 
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(2.18)  E[exp{it^Y^  +  it2Y2  )  ]  -  E[exp{ it^(I^Z^+»rj^Xj^^+it2 (l2Z2+^r2X2) )  ] 


(2.19)  E[exp{it^G^  +  it2G2)]  -  E[exp{ it^(I^K^+M^)+it2(l2K2+M2) ) ] 


where 


i  -  and  -®  <  t^, 


t2  <  «. 


Proof:  We  prove  the  lemma  for  the  exponential  case,  the  geometric  case  being 
similar.  Write 

Ni  N2  N2 

(I  X  I  X  )  -  (X  X  )+(x(N  >1)  I  X  .x(N  >1)  I  X,,) 

j-1  J  j-1  J  ^  j-2  J  j-2 

St 

where  x  denotes  the  indicator  function.  Clearly  > ’’’2^2^ 


we  are  left  to  show 


that  (x(N.,>l)  5]  X.,.,  x(N,,>l)  X  X.^.)  has  the  same  distri- 

^  j-2  2  j_2 


bution  as  , I2Z2)  .  Note  that 

Ni  N2 

(2.20)  P{x(N^>l)  I  X  <  X  x(N,>l)  I  X_.  <  X-} 

j-2  ^  ^  j-2 


1  1  "1  ^'2 

-I  I  P{X(N,>1)  X  X  <  X  .  x(N,>l)  X  X_,  <x,, 
i  -0  i  -0  j-2  j-2  ^ 


X(N^>1)  -  x(N2>1)  -  i2l 


Now,  for  (i^,i2)  -  (1,1)  in  (2.20)  we  have 


Ni  N2 

(2.21)  P{^X^  X^j  <  x^,  X^  X2J  <  X2.  Nj^>l,  N2>1) 
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00  00  1  “2 

-I  I  P(  I  X  <x  I  X  <  X  )P{N  -n  N,-n„)  . 

nj^-2  j-2  ^  ^  j_2  ^  ^  112  2 

But  P{N^-n^,  N2-n2}  -  P{N^-n^-l,  N2-n2-l)  P{Nj^>l,  N2>1),  so  that  (2.21)  is 
equal  to 


•i.  £. 

P{  X  X  <x  ,  X  X,,  <  X-)P{I  -1.  I-l) 
j-1  ^  J-1  ^  1  2 


-  P{Z^  <  x^,  Z2  <  X2)P{I^-1,  I2-I) 


For  (2.20)  we  have 


,22)  P{  X^  X^j  <  x^,  0  <  X2.  N^>1,  N2-I) 


00  n 


-  I  P{  I  X.  <x.,)P{N.-n,N,-l). 
n-2  j-2  J-J  1  1  2 


It  is  easy  to  check  via  (2.3),  that  P{N^-n.  N2-I)  -  Pio^Poo'^Pq!^  ^PlO'^Pll^’^'^ 
-  P(N^-n-l)P{N^>l,N2”l) ,  so  that  (2.22)  is  equal  to 


P{  <  x^}P{Ij^  -  1,I2  - 


-  P{Z^  <  x^)P{I^  -  l.l^  -  0). 


The  case  when  (i^,i2)  -  (0,1)  and  (ij^,i2)  “  (0.0)  in  (2.20)  follow  similarly. 


3.  The  General  BEAR(m)  and  BGAR(m)  Models 

In  this  section  we  construct  two  classes  of  AR  sequences  of  bivariate 
random  vectors.  In  each  class  the  sequences  are  labeled  by  the  parameter  m. 

We  denote  the  first  class  of  sequences  by  {X(m,n)  -  (X^(m,n) ,X2(ni,n) ) ' , 
n  -  0,1,...)  and  the  second  class  of  sequences  by  {G(m,n) 

-  (Gj^(m,n)  ,G2(ni,n) ) '  ,n  -  0 , 1 ,  .  .  .  )  ,  m  -  1 , 2 ,  .  .  .  .  We  show  that  the  random  vec¬ 
tors  X(m,n)  and  G(m,n)  have  bivariate  exponential  and  geometric  distributions, 
respectively,  with  mean  vectors  that  do  not  depend  on  m  or  n.  Then  we  discuss 
the  association  property  for  any  finite  number  of  random  variables  belonging 
to  one  of  the  two  AR  classes.  We  conclude  this  section  by  discussing  the 
autocorrelation  structure  of  the  two  classes  of  sequences.  Throughout,  n 
ranges  over  the  nonnegative  integers  and  i  assumes  the  values  1  or  2 .  For 
notation  simplicity  we  suppress  the  parameter  m  since  it  is  fixed  throughout 
the  section. 

First  we  construct  the  class  of  BEAR(m)  sequences.  Some  notation  is 
needed. 


Notation  3.^.  Let  ®  , jr2 (n)  ®  (0,1),  and  let  B(n)  be  a 

2x2  diagonal  matrix  with  B(n)  -  diag{jrj^(n)  ,»r2(n)  )  .  Further  let 
E' (n)  -  (Ej^(n) , E2(n) )  be  a  sequence  of  independent  bivariate  exponential  ran¬ 
dom  vectors  with  mean  vector  (Aj^^,A2^)',  let  e^  be  a  m-dimensional  vector  with 
component  j  equal  to  1  and  the  other  components  equal  to  0,j  -  l,...,m,  and 
let  0  denote  the  m-dimensional  zero  vector.  Finally  let  I ' (n)  - 
(I^(n,  1)  ,  .  .  .  Ij^(n,m)  ,  l2(n,  1) ,  .  .  .  , l2(n,m)  )  be  a  sequence  of  2m-dimensional 
independent  random  vectors  with  components  assuming  values  1  or  0  independent 
of  all  E(n) ,  and  let  A(n,q)  be  a  2x2  random  diagonal  matrix  with  A(n,q)  - 


diag{Ij^(n,q)  ,l2(n,q) )  ,  q  -  1 . m. 


m 


^4*1  *  »  t 


lu 

k 

'•i 

*3 


I 

I 


^  iJ 

»*!S 
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We  assume  that  for  i  -  1,2, 


(3.2)  I  P{(I  (n,l) . I,(n,m))  -  e.)  -  1-w/n), 

j_l  ^  J  ■« 


and  that 


(3.3)  P{(I^(n,l) . I^(n,m))  -  O')  -  »r^(n) 


We  define  the  BEAR(m)  sequences  as  follows. 


(3.4)  X(n)  - 


E(n)  n  -  0, . . . ,m-l 

m 

X  A(n,q)X(n-q)  +  B(n)E(n)  n  -  m,m+l,... 
q-1 


Next  we  construct  the  class  of  BGAR(m)  sequences.  Some  notation  is  needed. 

Notation  3.5.  Let  ,p2 ,aj^(n)  ,a2(n)  6  (0,1)  such  that  p^  <  o^(n),  let  N' (n)  - 
(N^(n) ,N2(n) )  be  a  sequence  of  independent  bivariate  geometric  random  vectors 
with  mean  vectors  (a^(n)p^^,  a2(n)p2^),  respectively,  and  let 
M' (n)  -  (Mj^(n) ,M2(n) )  be  a  sequence  of  independent  bivariate  geometric  random 
vectors  with  a  common  mean  vector  (Pj^^,P2^)  independent  of  all  N(n) .  Further 

let  J' (n)  -  (Jj^(n,  1) ,  .  .  .  , Jj^(n,m)  ,  J2(n,l) . J2(n,m))  be  a  sequence  of  2m- 

dimensional  independent  random  vectors  with  components  assuming  the  values  1 
or  0,  independent  of  all  N(n)  and  M(n) .  Let  C(n,q)  be  a  2x2  random  diagonal 
matrix  with  C(n,q)  -  diag{  Jj^(n,q)  ,  J2(n.q)  ) ,  q  -  1 . m. 

We  assume  that  for  i  -  1,2, 


(3.6) 


I  P{(J.(n,l) . J»(n,m))  -  e, )  -  l-a-(n) 
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and  that 


(3.7)  P{J/n,l) . J.(n.m))  -  O')  -  a  An) 


We  define  the  BGAR(m)  sequences  as  follows. 


(3.8)  G(n)  - 


M(n)  n  -  0,  .  .  .  ,m-l 

m 

5]  C(n,q)G(n-q)  +  N(n)  n  -  m.m+l,.... 
q-1 


Next,  we  show  that  X(n)  and  G(n)  have  bivariate  exponential  and  geometric 
distributions,  respectively. 


Lemma  3.9.  For  n  -  0,1,...,  X(n)  [G(n)]  has  a  bivariate  exponential 
[geometric]  distribution  with  mean  vector  (Aj^^,A2^)  [  (Pj^^  ,P2^)  ]  • 

Proof:  We  prove  the  results  of  the  lemma  by  an  induction  argvunent  on  n.  For  n 
-  0,...,m-l,  the  results  of  the  lemma  follow  by  the  definition  of  X(n)  and  of 
G(n) .  Let  us  assume  that  the  results  of  the  lemma  hold  for  all  nonnegative 
integers  that  are  less  than  or  equal  to  r,  r  >  m-1,  and  prove  that  the  results 
of  the  lemma  hold  for  r+1. 

t 

Let  E'  -  (£^^,£2)  [M  -  (M^,M2)]  be  a  bivariate  exponential  [geometric]  random 
vector  with  mean  vector  (Aj^^,A2^)  [(Pj^^,P2^)]  independent  of  all  E(n) 

[N(n)  and  M(n) ] .  Then,  by  the  induction  assumption,  we  have  for  i  -  1,2  that 


st  •  ''  P-  l-*'^(r+l) 

“  ■  x^(r+l)E^(r+l)  ,  w.p.  )r^(r+l) 

and  that 


V.W. 


+  N^(r+1).  w.p. 


l-a^(r+l) 


St 

G^(r+1)  - 


N^(r+1) 


,w.p.  a^(r+l), 


where  w.p.  stands  for  'with  probability'.  It  is  easy  to  check  that  X^(r+1) 
[G^(r+1)]  has  an  exponential  [geometric]  distribution  with  mean  [p^^] • 
Consequently  the  results  of  the  lemma  follow.  | | 


Now  we  consider  the  association  of  any  finite  collection  of  the  X^(n)'s 
and  of  the  G^(n)'s. 


Lemma  3.  W.  Let  us  assume  that  for  j 
X2(j)  in  (3.4)  are  associated;  let  n^^ 
and  let  6  (1,2),  r  -  1,2,.. 

-•  1 . r,  are  associated. 


“  0 . m-1,  the  random  variables  X^(j), 

<  n2  <. . .<  n^,  be  nonnegative  integers 

Then  the  random  variables  X.  (n  ) ,  q 

i  q  ^ 

q 


Proof;  Let  T2j  -  X2(j-1)  and  let  T2j.i  ”  j  ■*  1.2 .  To  prove  the 

result  of  the  lemma  it  suffices,  by  Barlow- Proschan  (1981,  Pj^,  p.  30),  to  show 
that 


(3.11)  the  random  variables  Tj^,...,T^  are  associated  for  all  r-1,2 . 


We  prove  (3.11)  by  an  induction  argument  on  r.  For  r  <  2m  (3.11)  follows  by 
the  lemma  assumption  and  by  Barlow- Proschan  (1981,  P^,  p.  30).  Let  us  assume 
that  (3.11)  holds  for  r,  r  >  2m  and  prove  that  (3.11)  holds  for  r+1. 


non- 


By  (3.4)  the  conditional  random  variable  . T^  is  stochastically 

decreasing  in  T^,...,T^.  Thus  by  Barlow- Proschan  (1981,  Lemma  4.8,  p.  147), 
there  is  an  r+1  argument  function  h,  nondecreasing  in  each  argument,  and  a 
random  variable  U  Independent  of  Tj^,...,T^,  such  that  (T^^ . ^r+1^ 

St 

-  (T^,...,T^,  h(U,Tj^,  .  .  .  ,T^))  .  By  Barlow- Proschan  (1981,  P2 ,  p.  30),  U  is 


I.V 
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associated  and  hence  by  Barlow- Proschan  (1981,  P^,  p.  30)  the  random  variables 


U,T^ . are  associated.  Consequently  by  Barlow- Proschan  (1981,  P^ ,  p. 

30),  {T^ . ^r+1^  associated.  || 


Using  a  similar  proof  we  obtain  the  following. 


Lenona  3.1^.  Let  us  assume  that  for  j  —  0 . m-1,  the  random  variables  Gj^(j), 

G2(j)  in  (3.8)  are  associated.  Let  n^^  <  n2  <.  .  .<  n^,  and  r  - 


1,2,...,  be  as  in  Lemma  3.10.  Then  the  random  variables  G.  (n),  q-l,...,r 

i  q  ^  ’ 

q 


are  associated. 


Finally,  for  each  class  of  sequences,  we  compute  the  autocorrelation 


functions  in  the  case  when  the  marginal  processes  are  stationary. 


For  the  exponential  models,  put  *'j(n)  -  all  n;i  -  1,2,  and  let 


P(l,(n,q)  -  1)  -  .^,(q),  Jt-  1,2;  q  -  1 . m 


such  that 


(i)  0/q)>O,  and  (ii)  I  ^,(q)  -  1-n „  i.  -  1,2, 


as  specified  in  equation  (3.2).  Define 


(k)  -  Corr(X^(n) ,X^(n+k) ) ,  i-1,2;  n  -  m,m+l , . . . ;k-l , 2 , 


Then 


(3.13)  p^  (k)  -  ^^(l)pj^  (k-l)+0 j(2)p.j,  (k-2)+.  .  .+(^B(m)p.^  (k-m) 


V 


with  Var{X^(n))  -  a'^,  i-1,2. 


For  the  geometric  models,  put  c«^(n)  -  a^,  all  n;  i-  1,2,  and  suppose  that 


I 


IWWJWIMUUIW W  WWMLIWWIW WIW  W  V  ">*  ■JTWTW  W  FUTTUT 1*  'WW  W  ".TV  "V  rv  n# ’V»rj«  -wt-TP.-v^rwijvjv'JV'jv'jwT^ 
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P{J^(n,q)  -  1)  -  7^(q).  1.2:q  -  1 . m 

such  that 

m 

(i)  7^(q)>0,  and  (ii)  J]  7^(q)-l-a^,  i  -  1,2, 

q-1 

as  specified  in  equation  (3.6).  Define 

(k)  -  Corr{G^(n) ,G^(n+k) ) ,  i-1 , 2 ;n-m,m+l , . . . :k-l , 2 . 

Then 

(3.14)  pg^(k)  -  7^(l)PQ^(k-l)+7^(2)/>g^(k-2)+...+7^(m)Pg^(k-m) 
with  Var(G^(n))  -  (l-p^)p'^.  t  -1.2. 

Evidently,  the  marginal  correlation  structures  of  the  bivariate  exponen¬ 
tial  and  geometric  sequences,  as  given  in  (3.13)  and  (3.14),  respectively,  are 
similar  to  that  of  the  Gaussian  AR(m)  process.  We  note  that,  in  general,  even 
when  the  marginal  processes  are  stationary,  the  joint  process  is  not  station¬ 
ary.  This  is  easily  seen,  for  example,  by  letting  m-1  in  (3.4)  with 
x^(n)  -  all  n;i  -  1,2;  choosing  E(n)  to  be  an  iid  sequence  of  random  vec¬ 
tors  where  Ej^(n)  and  E2(n)  are  iid  exponential  random  variables  for  all  n,  and 
letting  I(n)  be  an  iid  sequence  of  random  vectors  for  which  P{Ij^(n)-l, l2(n)-l! 

-  p^^  (1-x^)  (l-x^)  .  A  simple  computation  shows  that  Cov{Xj^(l)  ,X2(1) ) 

Cov{Xj^(2)  ,X2(2)  )  in  this  example. 

In  the  next  section,  we  develop  models  in  which  the  joint  processes  are 


also  stationary. 


4.  The  Stationary  BEAR(l)  and  BGAR(l)  Models 

In  this  section  we  consider  special  cases  of  the  BEAR(m)  and  BGAR(m) 
models  given  in  Section  3  in  which  the  joint  processes  are  stationary. 
Throughout  this  section  we  put  m  -  1,  assume  that  and  do  ^ot  vary 

with  n,  and  put  more  structure  on  the  E(n) ,  M(n) ,  and  N(n)  sequences.  We  show 
that  for  these  models,  the  bivariate  distributions  of  X(n)  and  of  G(n)  have  a 
form  of  the  type  studied  by  Arnold  (1975).  By  selecting  the  E(n),  M(n) ,  and 
N(n)  sequences,  as  defined  in  Section  3,  appropriately,  we  can  obtain  well 
known  bivariate  distributions.  For  a  stationary  BEAR(l)  model  we  obtain  the 
(i)  Marshall-Olkin  (1967),  (ii)  Downton  (1970),  (iii)  Hawkes  (1972),  and  (iv) 
Paulson  (1973)  bivariate  exponential  distributions.  For  a  stationary  BGAR(l) 
model  we  obtain  the  (i)  Esary-Marshall  (1973),  (ii)  Hawkes  (1972),  and  (iii) 
Paulson-Uppuluri  (1972)  bivariate  geometric  distributions.  We  conclude  this 
section  by  computing  the  autocovariance  matrices  for  each  model. 

First  we  present  the  exponential  case.  Some  notation  and  assumptions  are 
needed. 

Let  m  -  1  and  let  us  assume  that  (Ij^(n,l),  l2(n,l)),  given  in  (3.1),  is  an 
i.i.d.  sequence  of  bivariate  random  vectors. 

For  simplicity  of  notation,  denote  by  ^  and  let 

Pij  -  P{(I^(n,l),l2(n,l))  -  (i,j)),  i.j  -  0,1.  Note  that  by  (3.2)  and  (3.3) 

(4.1)  P^Q  +  P^^  -  l-,r^.  Pq^  +  p^^  -  1-n^. 

Further  let  (Nj^,N2)  be  a  bivariate  geometric  random  vector  with  parameters 
Pij ,  i, j  -  0,1  given  by  (2.3)  ,  and  let  E(r) ,  r  -  ±  1,  ±  2,  .  .  .  ,  be  an  i.i.d. 
sequence  of  bivariate  exponential  random  vectors  with  mean  vector  (Aj^^,A2^) 
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independent  of  (Nj^,N2)  and  all  (Ij^(n,l),  l2{n,l)).  Note  that  by  Lemma  2.14 


N, 


N, 


1  *’2 

(  X  X  ®  bivariate  exponential  random  vector  with 

j-1  j-1 

mean  vector  (A^^,A2^).  We  assume  that 


N, 


N, 


1  2 
(4.2)  E(0)  -  (  X  X  E  (-j))'. 

j-1  j-1 

Define  A(n)  to  be  the  2x2  diagonal  random  matrix  A(n)  -  diag{ I^(n, 1) ,  l2(n,l) 
and  B  to  be  the  2x2  diagonal  matrix  B  -  diag{ , ^2 } .  The  stationary  BEAR(l) 
model  is  defined  as  follows. 


(4.3) 


X(n)  - 


fE(0)  n-0 

A(n)X(n-l)  +  BE(n)  n-1,2, 


We  now  state  and  prove  a  characterization  of  X(n) 


•J 

K 


Lemma  4.4  Let  X(n)  be  defined  by  (4.3).  Then  for  n  -  0,1,..., 


where  E(0)  is  given  in  (4.2). 


St 

X(n)  -  E(0) 


Proof :  We  prove  the  result  of  the  lemma  by  an  induction  argximent  on  n.  By 
definition  the  result  of  the  lemma  holds  for  n-0.  Let  us  assume  that  the 
result  of  the  lemma  holds  for  n,  n  >  0. 

Note  that 


Ni  N2 

E'(0)  -  (  X  X  ’T  E-(-j)) 

j-1  j-1 


Ni  ..2 

(X(N^>1)  E^(-j),  X'(N2>1) 


^2 

I 

j-2 
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where  x( ’ )  denotes  the  indicator  function,  and  that  the  two  summands  are 
independent  random  vectors . 


Now  by  Lemma  2.17,  (4.1)  and  the  induction  assumption 


N, 


N. 


’1  2 
(X(N.>1)  »r  I  E.  (-j),x(N,>l)  w  I  E  (-j)) 
^  J-  j_2  ^  j-2 


St 


N, 


N, 


1  "2 
-  (I.(n,l)ir  I  E  (-j),I  (n,l)ff,  I  E  (-j)) 
^  j-1  ^  ^  j-1  ^ 


St 


-  (I^(n,l)X^(n),l2(n,l)X2(n)). 


Further  by  the  definition  of  E(r) 


St 


(»r,E.,(-l),  ir  E-(-l))  -  {w.,E.,(n),  7r,E„(n)) 


Since  the  random  vectors  (Wj^E^(n),  jt2E2(’^))  (I^(n,  l)Xj^(n)  ,  l2(n,  l)X2(n) ) 

are  independent  we  have  by  (3.4)  that 


*^1 

I 

j-1 


'1  ^^2  St 


+  (I^(n,l)Xj^(n).  l2(n,l)X2(n))  -  X' (n+1) 


Consequently  the  result  of  the  lemma  follows.  | 


In  the  following  remark  we  show  that  many  interesting  bivariate  distributions 
are  possible  in  Lemma  4.4.  See  Block  (1977)  for  further  details. 


■ 


a 


Remark  4.5.  (a)  Let  E  be  an  exponential  random  variable  with  mean  9, 

0  <  9  <  (A^+A2)  '^,  let  ir^  -  A^^ ,  and  let  E(l)  -  ^2^^)  ■ 

(i)  If  Pqo  “  Pqi  “  ’^l’  PlO  “  ^2  Pll  “  ^  resulting 

X(n)  has  independent  components. 

(ii)  Let  bj^,b2,b2  be  nonnegative  real  numbers  such  that  Ai  -  bi  +  b^2  and 

Aj  -  +  bjj.  If  -  «b^2.  -  »b^  p„j  -  «bj  and  p^^  -  1- Xb^+b^+b^^, 

then  the  resulting  X(n)  has  a  Marshal 1-Olkin  (1967)  bivariate  exponential  dis¬ 
tribution. 

(b)  Let  E(l)  have  independent  components: 

(i)  If  pj^j^  -  7(1+7)'^.  0  <  7,  p^Q  -  Pqj^  -  0,  and  p^^  -  (l-t-7)*^  then  the 
resulting  X(n)  has  the  Downton  (1970)  bivariate  exponential  distribution. 

(ii)  With  different  requirements  on  p^^^  i,j  -  0,1,  the  resulting  X(n)  has  the 
Hawkes  (1972)  and  the  Paulson  (1973)  bivariate  exponential  distribution.  See 
Block  (1977)  for  details. 

Next  we  present  the  geometric  case.  Some  notation  and  assumptions  are 
needed. 

Let  m  -  1  and  let  us  assume  that  (J^(n, 1) , J2(n, 1) ) ,  given  in  (3.5),  is  an 
i.i.d.  sequence  of  bivariate  random  vectors. 

For  simplicity  of  notation  denote  by  a^,  i  -  1,2,  and  let  p^^j  - 

P((Jj^(n.l).J2(n,l))  -  (i.j)},  i.j  -  0,1.  Note  that  by  (3.6)  and  (3.7) 

(4.6)  Pio  +  I’ll  ■  ^-“r  Pqi  *  Pii  ■ 

Further  let  (Q^,Q2)  be  a  bivariate  geometric  random  vector  with  parameters 
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Pj^.,  i,j  -  0,1,  given  by  (2.3).  Let  Pj^.P2  6  (0,1)  such  that  p^  <  ,  and  let 


N(r) ,  r  -  ±1,±2,...,  be  an  i.i.d.  sequence  of  bivariate  geometric  random  vec¬ 


tors  with  mean  vector  (a^p^^ ,Q!2P2^)  independent  of  (0^,0^)  and  all  (J^(n,l) 


Ql  Q2 

J.,(n,l)).  Note  that  by  Lemma  2 . 14  (  X  N,  (-j),  X  N-(-j))  is  a  bivariate 
2  j-1  j-1 

geometric  random  vector  with  mean  vector  (p^^  .p^  ).  We  assume  that 


^1  ^2 

(4.7)  M(0)  -  (  X  N  (-j),  X  N  (-j))' 

j-1  j-1 


Define  C(n)  to  be  the  2x2  diagonal  random  matrix  C(n)  -  diag( Jj^(n, 1) , J2(n, 1) ) 


The  stationary  BGAR(l)  model  is  defined  as  follows. 


(4.8)  G(n)  - 


C(n)G(n-l)  +  N(n)  n-1,2, . . . 


We  now  state  and  prove  a  characterization  of  G(n) . 


Lemma  4.9.  Let  G(n)  be  defined  by  (4.8).  Then  for  n  -  0,1,..., 


G(n)  -  M(0) 


where  M(0)  is  given  in  (4.7) 


Proof :  We  prove  the  result  of  this  lemma  by  using  a  similar  argument  to  the 


one  used  in  the  proof  of  Lemma  4.4.  || 


In  the  following  remark  we  show  that  many  interesting  bivariate  distributions 


are  possible  in  Lemma  4.9.  See  Block  (1977)  for  further  details. 


Remark  4.W.  (a)  Let  ”  Pi  **2  ”  ^2  N(l)  -  (1,1): 


(i)  Let  ^i'^2’^12  ^  1®*-  I'**!  “  ^1^12’  ^'“2  “  ^2^12'  l®t 


I.fl  <«•  *>l 


1*4.  a  4.tf 'A.A*%aA  4&l''i_A'4>4 


Pll  -  PlO  -  'l'l-''2''l2'  Poi  ■  <1-*1>'’2*12  “'*  Poo  -  ^-Pll-Poi-Pl 


Then  the  resulting  G(n)  has  the  Esary-Marshall  (1973)  bivariate  geometric  dis¬ 


tribution  in  the  narrow  sense. 


(li)  If  ^2.1  *  (f  (f  "*^2^  *  ^01  "  *  ^10  "  (f  ^00  **  ^1^2 


a,Q!„  then 


the  resulting  G(n)  has  independent  components. 


(iii)  With  no  special  requirements  on  the  p_  i,j  -  0,1,  the  resulting  G(n) 


has  the  Esary-Marshall  (1973)  bivariate  geometric  distribution  in  the  wide 


sense,  or  equivalently,  the  bivariate  geometric  distribution  due  to  Hawkes 


(1972)  . 


(b)  If  N(l)  has  independent  components,  the  resulting  G(n)  has  the  Paulson  and 


Uppuluri  (1972)  distribution. 


Finally,  we  give  the  autocovariance  matrices  for  the  ."stationary  BEAR(l) 


and  BGAR(l)  models.  Let  ”  Var{X(n))  be  the  variance  -  covar  ance  matrix  of 


X(n) ,  and  Y  -  Var{G(n)}  be  the  variance -covariance  matrix  of  G(n).  Note  that 

Cj  ““ 

Y  and  Y  are  independent  of  n  by  Lemmas  4.4  and  4.9,  respectively.  Define 

A  V 


r  (k)  -  Cov{X(n+k),X(n)),  k  -  0,1,2 .  and  F  (k)  -  Cov(G(n+k) ,G(n) ) ,  k 


0,1,2 .  and  note  that  rjj(O)  -  Yx  view  of  (4.3)  and 


(4.8),  it  is  easy  to  see  that  r.y.(k)  -  Arv(k-l)  and  F  (k)  -  CF_(k-l),  k  - 

A  A  (j  Cj 


1,2,...,  respectively,* where  A  and  C  are  the  2x2  diagonal  matrices  defined  by 


A  -  diag{ l-TTj^ ,  1-^2  )  and  C  -  diag{ l-a^^ ,  1-02 ) .  Hence,  for  the  stationary 


BEAR(l)  model  we  have 


(4.11)  Fj^(k)  -  A^Yy!  -r^(k).  k-0,1,2,. 


and  for  the  stationary  BGAR(l)  model  we  have 


(4.12)  Fg(k)  -  c'^Yg:  rG(-k)  -  Fg(k),  k-0,1,2, 


i’#>V*SV*4' 
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5.  The  BEARMA(ni^  ,ni^)  and  BGAR^(ra^  Models 

Using  the  results  of  Section  3  and  the  results  of  Langberg-Stoffer  (1985) 
for  MA  sequences,  we  construct  four  classes  of  ARMA  sequences  of  bivariate 
random  vectors.  In  each  class  the  sequences  are  labeled  by  the  parameters 
m^.m^.  We  denote  the  first  two  classes  of  sequences  by  {Z' (j ,m^ ,m2 ,n)  - 
(Z^(j ,m^ ,m2 ,n) , Z2(j ,m^ ,m2 ,n) ) ,  n  -  0,1,...),  j  -  1,2,  and  the  other  two 
classes  of  sequences  by  {L' (j ,m^,m2,n)  -  (L^( j , m^ ,m2 , n) , L2 ( j ,m^ ,m2 ,n) ) , 
n-  0,1,...},  j-  1,2.  We  show  that  the  random  vector  Z' (j ,m^ ,m2 ,n) 

[L' ( j ,m^ ,m2 ,n) ]  has  a  bivariate  exponential  [geometric]  distribution  with  a 
mean  vector  that  does  not  depend  on  j,mj^,m2,  or  n.  Then  we  discuss  the  asso¬ 
ciation  property  of  any  finite  number  of  random  variables  belonging  to  one  of 
the  four  ARMA  classes.  For  notational  simplicity  we  suppress  the  parameters 
m^  and  m2  since  they  are  fixed  throughout  this  section. 

First  we  construct  the  two  classes  of  BEARMA(mj^  ,m2)  sequences.  Some 
notation  is  needed. 

Notation  ^.1.  Let  X(n)  be  a  BEAR(mj^)  sequence  given  by  (3.4),  and  let  Y(n)  be 
an  m2-dependent  BEMA  sequence  as  given  by  Langberg-Stoffer  (1985,  p.  7), 
independent  of  the  X(n)  sequence.  Further  let  V' (n)  -  (V^(n) ,V2(n) )  be  a 
sequence  of  independent  bivariate  random  vectors  with  components  assuming  the 
values  1  or  0,  independent  of  the  X(n)  and  Y(n)  sequences  and  let 
P{V^(n)  -  1)  -  7r^(n),  0  <  7r^(n)  <  1,  f  -  1,2. 

We  define  the  two  BEARMA(mj^ ,m2)  sequences  as  follows. 

(5.2)  (Z^(l,n),Z2(l,n))  -  ( (l-,r^(n)  )Y^(n)  ,  (l-;r2(n)  )Y2(n) ) 


+  (Vj^(n)Xj^(n),V2(n)X2(n)) 


(5.3)  (Zj^(2,n),Z2(2,n))  -  ( (l-fl-2  (n)  )X2(n) ,  (n)  )X2  (n) ) 

+  (V^(n)Y^(n),V2(n)Y2(n)). 

Now  we  construct  the  two  classes  of  BGARMA(m^,m2)  sequences.  Some  notation  is 
needed. 

Notation  5.4.  Let  6  (0,1)  such  that  <5^  <  ^  ~  ^.2.  and  let 

G(l,n)  [G(2,n)]  be  a  BGAR(m^)  sequence  with  mean  vector  (^^^ , ’^2^2^^ 

[(^1  ’52^)]  given  by  (3.8).  Further  let  H(l,n)  [H(2,n)]  be  an  m2-dependent 
BGMA  sequence  with  mean  vector  (S^^,S2^)  [  (^^^5  ^  given  by  Langberg- 

Stoffer  (1985,  p.  8),  independent  of  all  G(l,n)  [G(2,n)].  Finally  let  y(n)  - 
(Uj^(n)  ,U2(n) )  be  an  i.i.d.  sequence  of  bivariate  random  vectors  with  com¬ 
ponents  assuming  the  values  1  or  0  independent  of  all  the  previous  random  vec¬ 
tors  such  that  P{U^(n)  -  1)  -  1-/9^,  -  1,2. 

We  define  the  two  BGARMA(mj^ ,m2)  sequences  as  follows. 

(5.5)  (L^(l,n),L2(l.n))  -  (G^(l,n) ,G2(l,n) 

+  (U^(n)H^(l,n),U2(n)H2(l,n)), 

(5.6)  (L^(2,n),L2(2,n))  -  (Hj^(2  ,n)  ,H2(2  ,n) ) 

+  (U^(n)G^(2.n),U2(n)G2(2,n)). 

Next  we  show  that  Z(j ,n)  and  L(j ,n)  have  bivariate  exponential  and 
geometric  distributions,  respectively. 

Lemma  5.7.  For  j  -  1,2,  and  n  -  0,1,...,  Z(j,n)  [L(j,n)]  has  a  bivariate 
exponential  [geometric]  distribution  with  mean  vector  (A’^,a‘^)  [(5.’^, 5'^)]. 


Proof :  By  Lemma  3.9  X(n)  [G(j,n)]  has  a  bivariate  exponential  [geometric]  dis¬ 
tribution.  By  Langberg-Stoffer  (1985,  Corollary  3.8),  Y(n)  [H(j,n)]  has  a 
bivariate  exponential  [geometric]  distribution.  Consequently  the  result  of 
the  lemma  follow  by  the  four  definitions  and  by  Lemma  2.10.  j  j 

Now  we  consider  the  association  property  of  any  finite  number  of  random  vari¬ 
ables  belonging  to  one  of  the  four  ARMA  classes.  We  assume  that  the  assump¬ 
tions  of  Langberg-Stoffer  (1985),  Lemmas  3.12  and  3.13)  are  satisfied.  We 
need  the  following  lemma. 

Lemma  5.8.  Let  S^, . . . ,S^,  T^, . . . ,T^  be  nonnegative  random  variables.  Let  us 
assume  that  . and  T^,...,T^  are  associated,  and  that  the  random  vec¬ 

tors  (Sj^ , . . . , S^)and  (T^,...,T^)  are  independent.  Then  the  random  variables 

StTi,...,S  T  are  associated. 

11  r  r 

Proof:  Let  T  -  (T^,...,T^),  let  W  -  (Sj^Tj^ ,  .  .  .  ,  S^T^) ,  and  let  f,g  be  two  nonne¬ 
gative  functions  each  with  r  arguments,  nondecreasing  in  each  argument. 

The  components  of  the  conditional  random  vector  W|T  are  nondecreasing  func¬ 
tions  of  the  associated  random  variables  Sj^,...,S^.  Thus,  by  Barlow- Proschan 
(1981,  Pj ,  p.  30),  the  components  of  W|T  are  associated.  Hence, 

E{cov(f(W),g(W))|T)  >  0. 

E{f(W)|T)  and  E{g(W)|T)  are  two  nondecreasing  functions  of  the  associated  ran¬ 
dom  variables  Tj^,...,T^.  Thus,  by  Barlow -Proschan  (1981,  P^,  p.  30),  the  two 
random  variables  E{f(W)|T)  and  E{g(W)|T)  are  associated.  Hence 

cov(E{f(W)|T),  E{g(W)|T))  >  0. 


Note  that 
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cov(f(W) ,g(W))  -  E{cov(f(W) .g(W)) |T) 

+  cov(E{f(W) |T>  E(g(W)|T)). 
Consequently  the  result  of  the  lenuna  follows.  || 


Lemma  5.9.  Let  us  assume  that  for  n  -  0,1,...,  V^(n)  and  V2(n)  are  associ¬ 
ated.  Let  n^  <  n_  <. . .<  n  ,  and  i- . i  ,  r  -  1,2 .  be  as  in  Lemma  3.10. 

12  r  1  r 

Then  the  (l,n  )  [Z^  (2,n  )],  q  -  1 . r  are  associated. 

q  q 

Proof :  By  Lemma  3.10,  X.  (n  ),  q  -  l,...,r  are  associated.  By  Langberg- 

q 


Stoffer  (1985,  Lemma  3.12),  the  random  variables  (n  ) ,  q  -  1, 


,r  are 


associated.  By  our  assumption  and  by  Barlow- Proschan  (1981,  P^,  p.  30), 
V  (n  ) ,  q  -  l,...,r  are  associated. 

q 


Thus,  by  Lemma  4.8,  (n  )  (n  ),  q  -  l,...,r,  and  (n  )  Y.  (n  ) ,  q  - 

1 . r,  are  associated.  Now,  Barlow- Proschan  (1981,  P^,  p.  30),  the  random 

variables  {X^  random  variables 

{Y^  (^q)’  ^2  associated.  Consequently  the  results 

of  the  lemma  follow  by  (5.2),  (5.3)  and  by  Barlow- Proschan  (1981,  P^ ,  p.  30). 

I  I 


Using  a  similar  argument  we  obtain  the  following. 

Lemma  5.W.  Let  us  assume  that  Uj^(O),  U2(0)  are  associated;  let 

Ui  <  n2  <. . .<  n^,  and  r  -  1,2 .  be  as  in  Lemma  3.10.  Then  the 

L,  (l.n  )  [L.  (2,n  q  -  1 . r,  are  associated. 

■'q  ^  "^q  '1''  J 

Langberg- Stoffer  (1985,  Section  4)  present  inequalities  and  probability 


bounds  for  the  bivariate  point  processes  related  to  the  bivariate  exponential 
or  geometric  MA  sequences.  We  note  that  all  the  results  given  by  Langberg- 
Stoffer  (1985,  Section  4)  hold  for  the  bivariate  point  processes  related  to 
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